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ABSTRACT 


NASA  has  recently  published  data  on  the  effects  of  tem- 
perature and  pressure  distortion  on  the  stability  of  axial- 
flow  compressors.   Theoretical  predictions  of  compressor 
instability  based  on  vorticity  do  not  explain  the  observed 
compressor  stability  of  the  NASA  tests.   Description  of 
distorted  flow  by  vorticity  is  extended  by  a  simplified  model 
which  incorporates  vortex  sheets  and  influence  coefficients. 
To  determine  the  upstream  influence  of  the  compressor  on 
flow  field  distortion,  a. model  using  parallel  compressor 
theory,  experimental  data  from  NASA  and  a  pressure  ratio 
equation  derived  by  Hill  and  Peterson  is  formulated.   The 
model  correctly  predicts  compressor  blade  loading  calculated 
by  Carta.   This  new  model  accounts  for  upstream  effects  of 
the  compressor  and  shows  that  the  vorticity  model  by  itself 
is  incomplete.   Incorporation  of  the  new  model  into  a  dis- 
tortion index  based  on  radial  vorticity  and  total  pressure 
change  is  proposed  and  discussed. 
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I.   INTRODUCTION 

In  applications  of  high-specific-energy  turbomachines , 
such  as  turbo-jet  engines  with  axial -flow  compressors,  it 
is  necessary  to  operate  the  compressor  as  close  to  maximum 
pressure  ratio  as  possible  without  inducing  unstable  per- 
formance known  as  surge.   A  phenomenon  of  axial  flow  com- 
pressors is  sensitivity  to  non-uniform  or  distorted  flow 
fields  when  operating  near  maximum  pressure  ratio.   As  air- 
craft are  designed  for  higher  speeds  and  more  comprehensive 
operating  envelopes,  the  possibility  of  flow  field  distor- 
tion at  the  compressor  face  is  increased.   At  supersonic 
speeds,  shock  effects  within  the  engine  intake  system  are 
significant.   At  high  altitudes  and  low  Mach  numbers  the 
effects  of  low  Reynolds  number  may  cause  portions  of  the 
flow  in  the  compressor  to  revert  to  laminar  flow.   Laminar 
flow  increases  the  tendency  toward  separation  which  accel- 
erates compressor  stall.   Aircraft  operations  at  high  angles 
of  attack  and/or  yaw  also  can  cause  distortions  of  the  flow 
field.   High  performance  military  aircraft  may  be  subjected 
to  all  of  the  above  distortions  in  addition  to  armament  re- 
lated temperature  distortions.   Forward  firing  armament 
could  also  cause  ingestion  of  foreign  gas  into  the  compressor 
and  is  discussed  in  Appendix  A.   Compressor  surge  is  not  re- 
stricted to  the  airborne  flight  envelope.   Ground  operations 
with  some  recently  developed  engines  have  resulted  in 
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destructive  surge  while  .taxiing.   Figure  1  shows  some  of  the 
general  areas  of  the  flight  envelope  where  large  stability 
margins  are  required.   Figure  2  shows  a  compressor  map  that 
illustrates  the  stability  margin.   The  associated  performance  - 
penalty  of  a  large  stability  margin  becomes  more  acute  as 
engine  performance  requirements  increase. 

This  situation  has  motivated  the  formulation  of  a  dis- 
tortion index  that  could  be  used  in  the  design,  development, 
and  component  matching  of  turbo-jet  engines  which  use  axial- 
flow  compressors.   Conceivably  such  an  index  could  be  used 
in  engine  control  systems  to  allow  operation  at  reduced 
stability  margins  without  the  danger  of  compressor  stall. 

Several  approaches  to  distortion  indices  have  been  tried. 
In  fact  there  has  been  a  proliferation  of  indices  based  on 
pressure  distortion  or  temperature  distortion  at  the  face 
of  the  compressor.   Most  data  collected  are  in  the  form  of 
total  pressure  or  total  temperature  measured  at  the  com- 
pressor face.   A  more  recent  approach  is  to  formulate  an 
index  based  on  vorticity.   Since  vorticity  can  be  derived 
from  either  temperature  or  pressure  data,  a  single  index 
appropriate  for  both  forms  of  distortion  is  possible. 

Farmer  [Ref.  1]  investigated  converting  total-pressure 
fluctuations  at  the  compressor  face  into  vorticity.   Iverson 
[Ref.  2]  proposed  a  method  of  transforming  total -  temperature 
contours  into  vorticity.   Shoemaker  [Ref.  3]  expanded  Farmer's 
theory  and,  using  a  computer  generated  map  of  vorticity,  de- 
rived an  index,  based  on  vorticity  alone,  tliat  predicted 
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Figure  1.   Flight  Envelopes  Showing  Areas  of  Potential 
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stall  in  selected  stall  events.   Research  conducted,  concur- 
rent with  this  paper,  by  Vann  [Ref.  4],  refines  the  Shoe- 
maker vorticity  index  and  applies  it  to  additional  stall 
events.   The  refined  vorticity  index  (Vann)  is  an  accurate 
stall  predictor  in  45  percent  of  the  cases  studied. 

Braithwaite,  Graber  and  Mehalic  [Ref.  5]  address  the 
effects  of  inlet  temperature  and  pressure  distortion  on 
turbo-jet  performance  and  provide  experimental  results  with 
which  the  vorticity  model  can  be  compared.   In  the  case  of 
total-pressure  distortion,  the  vorticity  model  is  a  good 
predictor  of  decreased  compressor  performance.   The  vorticity 
model  based  on  total -temperature  distortion  does  not  cor- 
relate well  in  all  cases  when  compared  with  Ref.  5. 

To  gain  further  insight  into  this  apparent  failure  of 
the  vorticity  model  to  correlate  with  experimental  results, 
this  paper  will  apply  the  vorticity  model  to  models  simu- 
lating the  conditions  of  Ref.  5.   In  addition,  an  extended 
model  that  includes  the  engine  effects  on  the  upstream  flow 
field  will  be  introduced  and  a  new  approach  to  a  distortion 
index  will  be  discussed. 
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II.   THE  NATURE  OF  THE  PROBLEM 


The  vorticity  approach  o£  Farmer,  Iverson  and  Fuhs 
[Ref.  6]  gives  the  following  relationships  for  pressure  and 
temperature  distortion: 


9P 


r   U'yP'r'   96* 
z  ' 


(1) 


,  ^  _M_   3T' 
'^r   2r'  dW 


(2) 


where  variables  with  superscripts  (primed)  are  non-dimen- 
sionalized  quantities  as  defined  below: 


r  =  r 


U  = 


w  =  to 


2 

U'a 

,  2a 
D 


P.  =  P'P. 


T   =  T  '  T 


D  =  compressor  diameter 


a  =  average  speed  of  sound 


w  =  fluid  vorticity 


P   =  average  stagnation  pressure 


T   =  average  stagnation  temper- 
ature . 


The  subscripts  6*  and  r  denote  circumferential  and  radial 
directions,  respectively.   A  complete  list  of  symbols  is  in- 
cluded on  pages  8  and  9. 

From  the  relationships  given  in  Ref.  6,  it  also  can  be 
shown  that  in  the  case  of  combined  pressure  and  temperature 
distortions 


U  0)   -  U  w 
z  r    1-  z 


RT   ^^t 


r    I^M^ 


rP    3  0* 


1  -^  ^  M' 


19T. 


P  r 
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am 


Order  of  magnitude  analysis  [Ref,  6]  indicates  that  U  cj   is 
of  negligible  size.   Reference  5  compares  the  effects  of 
pressure  and  temperature  distortion,  separately  and  in  com- 
bination, on  the  performance  of  J85-GE-13  turbo-jet  engines 
To  induce  total-pressure  distortions  on  the  face  of  the 
compressor,  screens  of  varying  solidity  and  180°  circumfer- 
ential extent  (semi-circle)  were  placed  one  diameter  upstre 
from  the  compressor.   The  total  temperature  of  the  air  was 
distorted  by  means  of  a  hydrogen  burner,  also  of  180°  extent, 
located  four  diameters  upstream  from  the  engine  inlet.   The 
burners  were  external  to  the  bellmouth.   The  test  procedure 
is  discussed  in  detail  in  Ref.  5.   Figure  3  shows  a  simpli- 
fied schematic  of  the  test  set-up  used.   Figure  4  is  a  com- 
pressor map  showing  the  degradation  of  performance  associated 
with  three  combinations  of  pressure  and  temperature  distor- 
tion at  the  engine  inlet.   From  Fig.  4  the  most  adverse 
effect  on  compressor  performance  is  observed  to  be  the  case 
of  temperature  and  pressure  distortions  super- imposed  (180° 
overlapped) .   Figure  5  shows  the  circumferential  profiles 
of  total  temperature  and  total  pressure  at  the  engine  inlet. 

Based  on  the  vorticity  approach  to  distorted  flow,  one 
would  expect  that  large  values  of  to  would  correlate  with 
large  loss  of  stall  margin.   Radial  vorticity  can  be  calcu- 
lated from  Eqn.  (3).   When  3P  /3G*  and  ST^/Se'-^  have  the  same 
algebraic  sign,  one  would  expect  the  most  unfavorable  case. 
Reference  to  Fig.  5(a)  shows  that  for  the  180°  overlapped 
case,  3P  /30*  and  3T  /36*  have  opposite  signs;  for  this  case 
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Figure  5.  Compressor  Inlet  Temperature  and  Pressure  Profiles 
witli  Combined  180°  Extent  Pressure  and  Temperature 
Distortions  (From  Ref .  5) . 
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(jj   is  small,  and  there  should  be  only  small  loss  of  stall 
margin.   For  the  0°  overlapped  case.  Fig.  5(c),  9P  /3e*  and 
8T  /36*  have  the  same  sign,  and  co   is  large.   The  0°  over- 
lapped case  should  result  in  large  loss  of  stall  margin. 
Observing  Fig.  4,  it  can  be  seen  that  the  opposite  result 
was  observed  experimentally. 
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III.   MODELS  OF  FLOW  WITH  PRESSURE  AND  TEMPERATURE  DISTORTION 

Effects  on  the  flow  field  at  the  face  of  the  compressor 
when  distortions  are  imposed  upstream  by  screens  (pressure) 
or  heating  (temperature)  will  be  examined.   Table  1  contains 
a  list  of  influence  coefficients  for  area  change  (dA/A) ,  heat 
addition  (dh  /c  T)  and  friction  (dCr/2)  that  will  be  used  in 
this  analysis.   These  influence  coefficients  are  adapted 
from  Shapiro  [Ref ,  7] . 

A.   PRESSURE  DISTORTION 

Figure  6  is  a  schematic  of  the  model  used  for  pressure 
distortion.   The  streamline  curvature  depicted  is  for 
graphic  purposes,  and  it  is  not  assumed  that  the  streamtubc 
change  actually  will  be  as  drawn.   The  screen  is  assumed  to 
be  of  180°  extent  (which  is  consistent  with  the  methods  of 
Ref.  5) ,  and  the  flow  at  station  1  is  assumed  to  be  uniform. 
From  continuity 


m.-|  =  m.^  =  p^U-A.   =  mass  flow  rate 
or  in  differential  form  after  logarithmic  differentiation 

a'  n         dp  ^  ^U  +  dA 


The  physical  geometry  requires  that 
lA.   ^  dA. 

a-'a  '  a-'b 


^1   =.^)  (4) 


The  assumed  uniform  static  pressure  distribution  requires 
that 
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AREA  CHANGE     HEAT  ADDITION       FRICTION 
dA  '^  ^t  ^   ^f 


A  .  ^p  T 


dU       1  1  Y  M^ 

U    M^  -  1         "   M2  -  1  r^M^ 


d£      Y  M^  Y  M^  [1  +  M'^  (y  -  1)]y  M^ 

p   ~  M2  -  1  M2  -  1  "        1  -  M2 


'^'^t  1 

L      0  - 

\  1  ^  ^-^  M 


2 


TABLE  1.   Influence  coefficients  for  area  change,  heat  addition 
and  friction. 


pU    p-'b  •  "■  ^ 


From  Table  1 


P-^B   M^  -  1    A-'B  "^  -* 


dP.   ^  -  Y  M'    dA.     n  +  M^(y-1)]y  M^   ^ 
P^A   M^  -  1     A-'a       1  -  M2  2   ■ 

(7) 

Combining  Equations  (4),  (5),  [6),  and  (7)  and  algebraic 
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manipulation  yields  the  following  relationship  for  area 
change : 


%A  =  [1  ^  M^  (Y  -  1)]  -~ 


dA. 
A''b 


(8) 


Using  similar  techniques  to  investigate  velocity,  it  can  be 
shovm  that 


dU 


dC, 


_1 dA.   .   Y  M^    ""f 

■  J  - 


U^A   M2  -  1    A^A    1  -  M2    2 


(9) 


and  combining  Equations  (8)  and  (9)  with  some  simplification 
reduces  to 

dC. 


dU. 
U^A 


[m^CI  +  y) -  1 
[   2(1  -  M2 


(10) 


From  Table  1 
dU 


U^B 


dA. 


1  -  M2    A^B 


(11) 


and  the  ratio  of  velocity  changes  can  be  expressed  as 


dU. 
dU. 

u-'b 


(Y  +  1)M'  -  1   ^ 
2       "  2 


dA 


1  -  M2     A'B 


) 


which,  using  Equation  (8),  reduces  to 


dU. 

u-'a 
u-'b 


(Y  ^  1)M^  -  1  ^  2.4  M^  -  1  ^  ^.j^.      .^2) 
1  +  (y  -  1)M2    1  +  0.4  M2    ^^^^  ^^^■' 


for  Y  =  1  •  4  . 

Equation  (8)  shows  that  for  all  real  valued  Mach  numbers, 
streamtube  A  will  increase  in  cross  sectional  area  (the  fric- 
tion coefficient  is  always  positive),  and  streamtube  B  must 
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decrease  a  like  amount  as  in  Equation  (4) .   For  subsonic 
flow,  Ug2  must  then  be  greater  than  Up, .   Equation  (12)  is 
finite  for  all  real  values  of  Mach  number  but  changes  sign 
when  M  =  0.6455.   For  values  of  M  less  than  0.6455,  Equa- 
tion (12)  is  negative  which  requires  the  velocity  change  in 
streamtube.  A  to  be  negative.   Conversely  when  M  >  0.6455, 
U.p  >  U, ,,   Quantitatively  M  may  be  less  than  0.6455  al- 
though transonic  compressors  have  axial  Mach  numbers  ap- 
proaching 0.65. 

B.   TOTAL  TEMPERATURE  DISTORTION 

Figure  7  is  a  schematic  of  the  total  temperature  distor- 
tion model.  The  illustrated  curvature  of  the  streamtube  is 
for  graphic  purposes  and  does  not  represent  an  assumed  pat- 
tern of  change.  At  station  2  the  static  pressure  distribu- 
tion is  assumed  to  be  uniform.  This  implies  that  station  2 
is  far  enough  ahead  of  the  compressor  so  that  there  is  no 
upstream  effect  of  the  engine. 

Using  Table  1 

P-'B    1  -  M2    A^B  ^      ■' 


,2   dh. 


dP^   ^   Y  M^    dA.   _   Y  M^   _t  f... 

P^A    1  -  M2    A-'a    1  -  M2  c  T   •  ^^■' 


Setting  Equation  (13)  equal  to  Equation  (14)  ,  because  of  the 
uniform  pressure  assumption,  then 

A^B    A^A      c  T   •  '■•'^^ 

P 
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Since 


dA.                dA. 

A-'b           a-'a  ' 

Equation    (15)    becomes 

dA.                dA^ 

A-'b           a-'a 

1  "Jl-t 

2CpT 

(16) 

With  the  addition  o£  heat,  -x)  u  decreases  and  — t-)a  increases; 

and  for  subsonic  flow,  U„-  must  be  greater  than  Un-,  . 

The  effect  of  heat  addition  on  the  velocity  is  investi- 
gated similarly 

dU    ^  _  _1 dA    ^  _1 ^ 

U-'A      1  -  M2   A^A    1  -  M2  c  T  ^^'■' 

P 

dU.   ^      1     dA.  ,^  J-. 

u-'b   "  1  -  M2  A-'b  •  ^^'^^ 

Incorporating  the  streamtube  area  changes  of  Equation  (16)  into 
Equations  (17)  and  (18)  ,  the  velocity  change  is  found  to  be 
the  same  for  both  streamtubes. 

U-'a    U-'b    2(1  -  M2)   c  T   •  ^^^■' 

C.   COMBINED  TEMPEMTURE  AND  PRESSURE  DISTORTION 

Figure  8  is  a  schematic  of  the  models  used  for  combined 
pressure  and  temperature  distortion.   Cases  will  be  considered 
for  zero  degrees  overlap  and  ]80°  overlap. 

From  Table  1  the  addition  of  heat  causes  a  loss  of  stag- 
nation pressure. 


dP. 

dh 

dh^ 

t 

-    Y   M    /2 

t 

Y   >4 

t 

P^ 

1      Y"l     .;i2 

c   T 

2 

c   C 

t 

1  +  ^fy-    M^ 

P 

P    t 
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In  the  case  of  0°  overlap 


dP. 


-  y  M^ 


^  ^A 


dC, 


(20) 


dP. 


) 


-  Y  M' 


dh. 


Pt  '^ 


P  t 


(21) 


The  net  area  change  of  streamtube  A  is,  using  Equations  (8) 
and  (16)  , 


J .  dC r    1   dh. 

A  ''aNET    ^^         lY  -LJ  H  J   ^     2   c  T 


(22) 


To  obtain  results  for  a  specific  case,  let  us  assume  that 
the  heat  and  screens  are  such  that 
dC,    ^  dh^ 


2     2  c  T^ 
P  t 


Consequently 


dP. 


dP. 


) 


) 


Pt  '^  \    '^ 


Equation  (22)  in  terms  of  T   and  introducing  dC^/2  in  terms 
of  heat  becomes 

'"l   dh^ 


dA  .      _  [1  +  (y-1)  mm 
A  -'ANET         2 


o   c  T. 
2    P  t^ 


1  + 


X^M^ 


1  dh. 


c  T, 
P  t 


dA  . 
A  ''aNET 


1  ^\ 


4  FTT 
p  t 


2    2 


(23) 


For  velocity  change  in  streamtube  A, 
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dU  .   ^    1    "  dA  .         Y  M^    ^ 
U  -"a    1  -  M^    A  ''ANET    1  -  M2    2   ' 

and  incorporating  Equation  (23)  simplifies  to 

dU  .   ^  1  +  2y  M'    ^  f 

U  -"b    (1  -  M2)  2    2    •  '■^^■^ 

In  streamtube  B, 

dU  .   ^   -1     dA  .         1    ^ 

U  ^B    1  -  M2    A  ^BNET    1  -  M^  c  T 

P 

dU  ,   _    1      dA  ,      ^1      '^^t 


IiJb      1-M2      aJ 


+ 


U  ^B    1  -  M2    A  ^ANET    1  -  M^   c  T 

P 

dU  .   _  5/2  +  (Y  -  1)M^   ^  .      . 

U  ^B        1  -  M2         2    •  ^   '' 

The  ratio  o£  velocity  changes  in  the  two  streamtubes  (for 
Y  =  1.4)  is 


)  V>  1   Ca.n;1\l2 


dU 

U  ^B    1.5  +  0.4  M'     rr^,^  roA  ^ 

dT—  ^  0.5  +  1.4-MF-  =  f (^^  f26a) 

u  -"a 


Examination  of  Equation  (26a)  shows  that  the  velocity  in 
streamtube  B  is  higher  for  all  subsonic  values  of  M.   Recall 
that  this  result  is  for  the  specific  case  of  equal  total 
pressure  in  streamtubes  A  and  B  which  resulted  from  the  ar- 
bitrary relation  of  dC^/Z  =  1/2'dh^/c  T^ . 
^  f      '     t   p  t 

For  the  case  of  180°  overlap  the  same  method  of  develop- 
ment is  used.   Figure  8(a)  is  the  model  used,  and  there  are 
three  stations  illustrated  for  the  purpose  of  adding  the 
solutions  of  each  distortion.   Assume  that  static  pressure 
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is  uniform  at  stations  2. and  3.   Adding  Equations  (19)  and 
(10)  and  again  considering  the  specific  case  of 


dC, 


1  ^\ 


2  2   c  T^ 

P  t 

the  velocity  change  in  streamtube  A  is 

dCr 


dU  . 

u  -'a 


4(1-M2) 


(1  +  2yM^) 


For  streamtube  B,  using  Table  1,  Equations  (8),  (16),  and 

(19)  the  velocity  change  is 

JTT        dC^ 

^)b  =  4TT^  [3  .  2M^(Y-1)] 

and  the  ratio  of  velocity  changes  is,  for  y   =  1.4, 


dU  . 

u  -"a 

dU  . 

u  -"b 


1  +  2  M' 


1  +  2.8M2 


3+2M2 (y-1)    3  +  0.8M2 


f(M) 


(26b) 


Examination  of  Equation  (26b)  shows  that  the  velocity  in 
streamtube  B  is  always  higher  than  the  velocity  of  stream- 
tube  A.   For  the  case  of  180°  overlap  of  distortions,  the 
total  pressure  will  not  remain  equal  in  both  streamtubes 
as  in  the  case  of  0°  overlap.   In  the  case  of  180°  overlap, 
the  change  in  total  pressure  for  the  distorted  streamtube 
(A)  will  be 


-2yM'    "^^f 
2 

for 

this 

specific    case    of 

dC^          1      ^^t 
2              2      CpT^ 

32 


Total  pressure  will  remain  constant  in  the  undistorted  stream- 
tube  (B)  . 

D.   HEAT  ADDITION  AT  SMALL  MACH  NUMBERS 

Using  the  temperature  distortion  model  of  Figure  9,  the 
effects  of  heat  addition  on  the  flow  field  will  be  made  more 
restrictive  by  limiting  the  entering  Mach  number  to  very 
small  (near  zero)  values.   This  is  accomplished  by  placing 
the  burner  outside  the  bellmouth,  where  velocities  are  small. 
This  analysis  also  models  the  tests  by  NASA  reported  in 
Ref.  5.   An  energy  balance  at  stations  1  and  2,  assuming 
adiabatic  conditions  between  these  stations,  gives 

^A^Va2  ^  l^iz)  ^  '^B^VbZ  ^  I"b2^   •   .    ^27) 

Note  that  station  1  is  immediately  downstream  from  the  burner. 
For  streamtubes  A  and  B  then 


and 


VbI  ^  I  "bI  =  Vb2  ^  ^  "b2  ^  VtB   •         .  ^29) 


Assume  uniform  P  at  stations  1  and  2,  m.  =  m„  and  isentropic 
flow  for  streamtubes  A  and  B  between  stations  1  and  2.   With 
the  previous  assumption  of  small  Macli  number,  conditions  at 
station  1  are  specified.   The  stagnation  pressure  for  both 
streamtubes  at  either  station  1  or  2  equals  the  ambient 
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Figure  9.      Temperature  nistortion  Model  for  Heat  Addition  at 
Low  Mach  Number. 
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pressure.   Static  pressure  for  both  streamtubes  is  equal  at 
station  1.   Due  to  heat  addition  T  .  >  T  g.   The  following 
isentropic  relationships  now  can  be  used: 

Y-1 
T..  =  T_  Ik^  1  (30) 


'A2    '^tA  IP 


tAl 


and 


T    =  T 
B2     tB  IP 


B2 


tBl, 


Y-1 
Y 


(31) 


Using  Equations  (30),  (31)  with  (28)  and  (29)  will  give  the 
relationship  of  velocity  for  streamtubes  A  and  B 


U^   =  2c  T 
^A2    ^  p^tA 


11^   =  2r  T 
^B2    ^  p  tB 


1  - 


Y-1 


Y-1 
\   Y 


Since 


K2 


tA 


T 


P_ 

P. 


) 


tB 


P_ 
P. 


)b  -  f(Y,M) 


(32) 


and  for  hydrogen  (fuel  for  burner  in  NASA  test  of  Ref.  5) 

oxidized  in  excess  air,  Ya  ~  Y^^.  then  M^2  ^  ^^B2'        ^^^^ 

Investigating  tlie  area  change  at  2  witli  the  same  assuinp 
tions 
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pAU)32  =  PAU)^2 


''B2    ^B2    A2 
and  using  Equation  (32) 


A 

p 

T 

/T 

A  2 

^B2 

^A2 

•^Hb 

^B2 

^2 

T 
^B2 

A,, 

Ab2     /r^ 


(34) 


Equation  (34)  is  restrictive  in  tliat  the  assumptions  limit  it 
to  the  case  o£  equal  mass  flow  in  streamtubes  A  and  B.   How- 
ever, parallel  compressor  theory  uses  a  mass -weighted  sum  of 
subcompressors  to  determine  the  operating  point  of  a  distorted 
compressor,  so  in  the  case  of  180°  sector  distortion  this  as- 
sumption is  not  unreasonable. 
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IV.  COMPARISON  OF  FLOW  MODELS  WITH  VORTICITY  THEORY  FOR  TEM- 
PERATURE DISTORTION 

The  analysis  of  Section  III(B),  total  temperature  dis- 
tortion, shows  that  the  velocities  of  the  two  streamtubes 
are  the  same,  thus  making  the  radial  vorticity  zero.   Set- 
ting to   in  Equation  (3)  equal  to  zero  and  solving  for  9P  /Se* 
yields  the  following  relationship: 


RT^ 

3P, 

90* 

-      S 

3Tj 

■     iji  M^  ■ 

'■Pt'^z 

Li  .  _Iji  H^J 

Simplifying  and  changing  T  to  T   gives 


RT^  dP^ 


Y-1 


c   dT^   ^ 
p    t   2 


M' 


dP^        T  dh^     c 
Pt       2   RT^  ^   Cp 


and  since  R/c   =(Y'iyY 


dP. 


1   c    dh, 

Y-1  _£_     t 

2     R   c  T 


P  t 


M' 


yM' 


dh. 


p  t 


(35) 


Equation  (35)  is  equivalent  to  the  influence  coefficient  of 
Table  1,  thus  verifying  vorticity  for  this  case. 

For  the  case  of  heat  addition  at  low  Mach  number,  a  more 
subjective  comparison  must  be  made.  Equation  (3)  simplified 
for  temperature  distortion  only  shows 
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c    3T^    I  ^  M^' 


All  values  in  the  right  hand  side  of  Equation  (36)  are  always 
positive  with  the  exception  o£  the  circumferential  total  tem- 
perature gradient.   For  the  velocity  profile  shown  in  Figure 
9,  w   is  positive  when  0*  =  180°  and  3T  /80*  is  also  posi- 
tive.  A  similar  comparison  at  6*  =  0°  also  shows  algebraic 
agreement  in  that  03   and  the  total  temperature  gradient  are 
both  negative. 
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V.   EXTENDED  MODEL  TO  ACCOUNT  FOR  COMPRESSOR  INFLUENCE 

A.  PARALLEL  COMPRESSOR  MODEL 

Reference  5  employs  the  parallel  compressor  model  as  a 
method  of  predicting  the  effects  of  distortion  on  stall 
pressure  ratios".   Figure  10  sho\^/s  a  schematic  of  the  parallel 
compressor  model  with  some  of  the  assumptions  used.   Figure 
11  shows  the  predicted  stall  points  of  the  compressor  for 
various  pressure  and  temperature  distortions  using  this 
model.   Figures  12  and  13  show  the  profiles  of  static  pres- 
sure, total  pressure  and  total  temperature  at  the  compressor 
inlet  and  outlets  for  both  pressure  and  temperature  distor- 
tion.  From  Figure  12  it  can  be  seen  that  static  pressure 
is  not  uniform  at  the  face  (inlet)  of  the  compressor.   This 
was  one  of  the  assumptions  used  in  the  development  of  Sec- 
tion III  of  this  paper.   The  development  of  an  extended 
model  in  this  section  will  consider  the  effects  of  non- 
uniform pressure  profiles. 

B.  PRESSURE  RATIO  EQUATION 

Hill  and  Peterson  [Ref.  8]  derive  a  compressor  pressure 
ratio  equation  that  predicts  the  total  pressure  ratio  (IT  ) 
across  an  axial  compressor  stage  as 

tan 


n 


32 


tan 


31 


(37) 


where  the  subscripts  1  and  2  refer  to  stations  before  and 
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FLOW 


STATION  2 


1.  SECTORS   C  AND   D   ARE  CONSIDERED  AS 

INDEPENDENT  PARALliL  COMPRESSORS 

2.  BOTH  SECTORS  FOLLOW  SAISf  ^^AP  AS 

UNDISTORTED  COMPRESSOR 

3.  NO  CROSSFLOW  BETWEEN  PARALLEL  COM 

PRESSORS 

4.  BOTH  EXIT  TO  A  COM/WN  STATIC  PRESS 


Figure  10.   Parallel  Compressor  Model  (From  Ref.  5) 
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Figure  11.   Application  of  Parallel  Compressor  Model  to  Predict 
Stall  Points  (I-rom  Rcf.  5). 
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after  the  compressor  stage.   Figure  14  shows  the  relationship 
of  the  angles  3,  and  B^  with  a  compressor  blade  row  at  wheel 
speed  Qr.      For  flow  in  unstalled  cascades  of  high  solidity 
(ratio  of  blade  chord  length  to  blade  spacing  in  the  cascade) , 
the  exit  angle  3-;,  is  relatively  constant  with  changes   of  3-,. 
Blade  loading  is  a  function  of,  or  can  be  correlated  to,  A3. 
If  Q>y    -  constant,  then  3^  has  a  significant  effect  on  the 
stall  characteristics  of  a  cascade.   Consider  Equation  (37) 
with  a  fixed  H   and  wheel  speed.   An  increase  in  temperature 
must  be  accompanied  by  an  increase  in  3-.  .   An  increase  in  3^ 
can  be  obtained  only  by  a  decrease  in  axial  velocity.   This 
is  shown  in  Figure  15.   Figure  11  shows  that  parallel  com- 
pressor theory  substantiates  the  decrease  in  axial  velocity 
for  this  case  since  corrected  weiglit  flow  is  proportional  to 
axial  velocity.   The  model  of  Section  III.D.  shows  a  veloc- 
ity profile  reversed  from  the  results  of  this  section;  see 
Figure  9.   An  assumption  of  Section  III.D.  was  uniform  static 
pressure  at  the  downstream  station.   The  model  in  fact  pre- 
cluded any  upstream  effect  of  the  compressor  itself  on  the 
flow  field. 

C.   PRESSURE  DISTORTION 

Figure  16  illustrates  the  profiles  of  static  and  total 
pressure?  static  and  total  temperature  and  axial  velocity 
as  derived  from  the  models  of  Section  III  and  as  measured 
at  the  face  of  the  compressor  during  the  NASA  tests  of  Rcf. 
5.   Flow  properties  at  the  constant  -  static -pressure  station 
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agree  with  those  illustrated  by  Figure  6,   According  to 
Figure  12,  the  static  pressure  increases  in  the  high-P 
streamtube  and  decreases  in  the  low-P   streamtube.   The 
flow  decelerates  in  the  high-P   streamtube  and  vice-versa 
for  low-P   streamtube.   Also  the  area  of  the  high-P   stream- 
tube  of  Figure  6  must  increase. 

The  screen  creates  a  vorticity  sheet  which  is  convected 
into  the  compressor.   This  sheet  separates  the  high-P   and 
low-P.  streamtubes.   As  predicted  by  Helmholtz's  vorticity 
laws,  the  vortex  sheet  is  a  flow  feature  which  persists  once 
it  is  created  [Ref,  9].   The  vortex  sheet  decays  by  means 
of  viscous  forces.   Figure  17  shows  a  model  that  uses  a 
"bent  vortex  sheet"  associated  with  the  change  of  velocity 
and  streamtube  area  for  pressure  distortion.   Specific  points 
(A,  B,  C,  D)  in  relation  to  this  vortex  sheet  are  indicated, 
and  the  change  in  axial  velocity  at  these  points  is  also 
shown  in  Figure  17. 

Due  to  bending  of  the  vortex  sheet,  the  flow  acquires  a 
Uq  component  in  the  direction  of  rotor  motion  at  points  A 
and  B  and  a  Uq  opposed  to  rotor  motion  at  points  C  and  D  of 
Figure  17.   This  has  an  influence  on  3-,  as  discussed  later. 

D.   TEMPERATURE  DISTORTION 

Figure  18  shows  the  flow  properties  at  three  stations: 
immediately  downstream  of  tlie  burner,  at  a  constant- static- 
pressure  location  within  the  duct  and  at  the  rotor  inlet 
station.   The  flow  properties  at  the  const  ant -static-pressure 
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station  within  the  duct  agree  with  Figure  9  of  Section  III 
and  figure  13 (a) . 

To  determine  the  upstream  influence  of  the  compressor, 
focus  attention  on  Figure  11(b)  which  is  the  parallel  com- 
pressor model  matching  Figure  18.   The  hot  streamtube, 
identified  as  1  in  Figure  11(b),  has  sharply  reduced  cor- 
rected weight  flow.   This  implies  low  axial  velocity  com- 
pared with  the  cold  streamtube,  number  2  in  Figure  11(b). 
Consequently,  the  velocity  profile  reverses  from  fast  in 
the  hot  streamtube  at  the  constant-static-pressure  station 
to  slow  in  the  hot  streamtube  at  the  rotor.   This  reversal 
is  attributable  to  upstream  influence  of  the  compressor. 
The  original  vorticity  sheet  generated  by  heat  addition  re- 
mains.  The  static  pressure  is  not  uniform  at  the  rotor, 
which  was  an  essential  element  of  the  models  of  Section  III. 
The  effect  of  the  compressor  must  be  to  change  the  areas  of 
the  streamtubes  by  an  increased  corrected  weight  flow  in 
the  undistorted  streamtube.   Figure  19  shows  the  bent  vor- 
tex sheet  as  applied  to  the  temperature  distortion  model 
and  the  velocity  changes  at  the  rotor. 

E.   EFFECTS  OF  THE  BENT  VORTEX  SHEET  ON  ROTOR  BLADE  INCIDENCE 
ANGLE 

Carta  and  Adamczyk  [Ref.  10]  postulate  that  it  is  not 

possible  to  prescribe  the  conditions  at  the  compressor  face 

a   priori    on  the  basis  of  upstream  distortion  profiles  alone. 

As  discussed  in  Ref.  ]0  tlic  upstream  conditions  are  strongly 

dependent  on  the  presence  of  the  blade  row  itself.   Carta  and 
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Adamczyk  present  a  graphic  comparison  of  angle  of  incidence 
distribution  as  it  differs  immediately  behind  the  distortion 
screen  (upstream)  from  that  at  the  blade  row  inlet.   This 
angle  of  incidence  distribution  of  Ref.  10  is  compared  in 
Fig.  20  with  a  semi -quantitative  analysis  of  inlet  flow 
angle  based  on  the  bent  vortex  sheet  model  for  pressure  dis- 
tortion.  This  analysis  is  discussed  in  Appendix  B. 

There  are  no  measured  or  predicted  angle  of  incidence 
distributions  available  for  comparison  v-zith  the  temperature 
distortion  model.   However,  the  encouraging  agreement  in  the 
case  of  pressure  distortion  prompts  a  similar  analysis  of 
temperature  distortion.   Figure  19  shows  the  change  in  axial 
velocity  at  the  points  A,  B,  B',  C,  D  and  E,  and  Figure  21 
shows  the  change  in  U„  at  these  points.   A  semi-quantitative 
analysis  similar  to  Appendix  B  gives  the  inlet  flow  angle 
(3-.)  distribution  shown  in  Figure  22,   The  area  above  the 
nominal  clean  inlet  value  of  B-,  contributes  to  stall.   The 
large  angles  of  3-.  in  Figure  22  lead  the  distorted  flow  sec- 
tor somewhat,  but  the  majority  of  the  high  blade  loading  due 
to  large  values  of  3-,  lies  in  the  distorted  sector.   Observ- 
ing Figure  11(b),  it  can  be  seen  that  the  compressor  with 
distorted  flow  stalls  first,  which  gives  qualitative  verifi- 
cation of  the  bent  vortex  model  for  temperature  distortion. 
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VI.   DISTORTION  INDEX  FORMULATION 

The  results  o£  Section  V  suggest  an  index  to  predict 
compressor  instability  or  stall  that  combines  total  pressure 
and  radial  vorticity.   Figure  21  indicates  that  the  area 
experiencing  the  highest  values  of  inlet  flow  angle  are  those 
that  also  would  experience  positive  radial  vorticity.   Since 
there  are  existing  data  and  digital  computer  programs  [Ref.  3] 
that  compute  radial  vorticity  as  well  as  AP  ,  then  it  is 
easy  to  search  the  compressor  face  for  these  values.   There 
are  several  possibilities  for  processing  these  values,  once 
they  are  calculated,  into  an  index.   One  method  would  be  to 
compare  the  sum  of  radial  vorticity  and  AP   for  a  specified 
sector  (percentage  of  compressor  face  area)  vv'ith  a  critical 
value.   Since  there  is  no  known  critical  value,  another  ap- 
proach that  eventually  might  lead  to  a  critical  value  is 
recommended.   In  Figure  23  the  total  pressures  across  the 
face  of  the  compressor  are  shown.   P.  is  the  average  total 
pressure.   If  P   -  P   is  positive,  then  the  effect  is  to 
increase  the  blade  loading  as  shown  in  Figure  20.   Positive 
radial  vorticity  also  increases  blade  loading  as  shown  in 
Figure  20(a),  and  as  in  Refs.  1-3  and  6.   The  most  signifi- 
cant effect  on  compressor  instability  or  stall  would  be  the 
largest  value  of  AP   and  radial  vorticity.   The  pressure 
change  can  be  normalized  with  P  ,  and  an  index  can  be  formu- 
lated by  summing,  over  a  sector,  the  combined  values  of  AP 

and  (jj  . 
r 
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k 
INDEX  =   E   (AP    +  u   ) 
i=l    ^i    ^i 


where  the  integer  k  represents  the  number  of  pressure  rakes 
to  be  searched.   By  rotating  or  stepping  the  sector  around 
the  array  of  probes,  it  is  possible  to  compute  the  maximum 
valued  sector  at  any  time.   By  plotting  the  maximum  value 
of  the  index  versus  time,  a  graph  similar  to  Figure  24 
could  be  analyzed  for  a  possible  critical  value.   A  time 
history  statistical  comparison  of  deviation  from  the  mean 
value  of  the  index  is  another  possible  test  to  be  applied 
and  is  done  by  Vann  [Ref.  4]  after  adapting  this  proposed 
index  to  existing  digital  computer  programs  and  comparing 
it  with  other  known  indices. 

It  is  implied  that  the  index  is  a  positive  number  at 
all  times.   There  is  nothing  in  the  formulation  that  re- 
quires the  index  to  be  positive,  and  if  the  sector  over 
which  the  index  is  computed  is  small  enough  then  it  may  in 
fact  be  negative.   Throughout  this  paper  the  term  stall  has 
been  applied  to  large  inlet  flov\?  angles  that  load  the  blade 
row  beyond  the  stall  limit  of  the  individual  blades.   If 
the  inlet  flow  angle  is  too  small,  a  phenomenon  known  as 
negative  stall  occurs.   Reference  8  discusses  the  plicnomenon 
of  negative  stall,  but  it  is  important  to  point  out  that  a 
negative  distortion  index  would  not  indicate  a  negative 
stall.   In  tlie  analysis  of  this  index  by  Vann  [Ref.  4], 
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a  sector  of  60°  was  used,  and  the  index  values  were  always 
positive.   Reference  4  addresses  the  size  of  sector  used  to 
compute  this  index  and  other  commonly  used  indices. 
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VII.   CONCLUSIONS 

Interpretation  of  distorted  inlet  flow  interaction  with 
a  compressor  based  on  vorticity  maps  [Refs,  1-3]  is  satis- 
factory as  long  as  the  assumption  of  constant  static  pres- 
sure is  valid.   Experiment  has  shown  that  constant  static 
pressure  is  not  the  case  at  the  compressor  face,  and  the 
vorticity  model  must  be  expanded  to  account  for  this.   The 
failure  of  the  vorticity  model  to  predict  the  conditions 
at  the  face  of  a  compressor  (and  thus  its  loss  of  stall  mar- 
gin) in  the  case  of  total  temperature  distortion  is  a  result 
of  not  considering  the  effects  of  the  compressor  on  the  up- 
stream flow  field.   The  upstream  influence  of  the  compressor 
causes  an  area  change  of  the  undistorted  streamtube  that  re- 
sults in  greater  velocity  changes  than  those  induced  by  the 
vortex  itself.   The  pressure  ratio  equation,  Equation  (37), 
and  parallel  compressor  theory  adequately  predict  trends  in 
compressor  stability. 

A  model  based  on  the  bending  of  the  vortex  sheet  gives 
results  comparable  to  parallel  compressor  theory  and  in  the 
case  of  total  pressure  distortion  qualitatively  duplicates 
Carta's  results. 

Carta's  stall  incidence  angle  ratio  suggests  a  distor- 
tion index  that  combines  radial  vorticity  and  total  pressure 
gradients  since  both  can  be  related  to  blade  loadings  that 
cause  stall.   Such  an  index  has  been  developed  and  lias  boon 
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tested  [Ref.  4]  against  other  distortion  indices  using  common 
data  for  comparison.   The  new  distortion  index  has  proven 
equal  to  or  better  than  the  distortion  indices  with  which  it 
was  compared  for  accuracy  in  predicting  compressor  stall. 
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APPENDIX  A 

Ingestion  of  hot  gas  from  forward  firing  armament  sys- 
tems has  more  than  temperature  effects  on  the  performance 
of  an  axial-flow  compressor.   The  characteristic  line  of 
corrected  engine  rotational  speed  (N//6  =  constant)  appears 
to  be  a  function  of  temperature  alone  since  6  is  defined 
as  the  ratio  of  actual  temperature  to  some  reference  tem- 
perature, usually  a  sea  level  standard  day  temperature. 
However,  the  reference  to  temperature  alone  is  a  convention 
that  has  evolved  from  the  fact  that  compressors  in  gas 
turbine  engines  for  aircraft  applications  use  only  air  as  a 
fluid  medium.   The  corrected  engine  rotational  speed  (N//0) , 
sometimes  identified  as  referred  speed,  is  a  similarity 
parameter  designed  to  ensure  comparability  of  data  taken 
under  varying  conditions  of  testing.   In  the  case  of  a  com- 
pressor blade  row  it  is  desirable  to  preserve  the  blade 
angle  of  attack.   Consider  two  "different"  compressors  as 
in  Figure  A-1.   If  the  angle  of  attack  is  to  remain  the  same 
for  both  compressors,  tlien  3  must  be  equal  to  3'.   6  is  the 
inlet  flow  angle  but  differs  from  tlie  angle  of  attack  only 
by  a  constant  related  to  blade  geometry  as  discussed  in  Ap- 
pendix B.   For  similarity 

„    Qr   f2'r' 
tan  3  ~ 


U   "  U 


and   M   must    equal    M'.       Tlicrefore 
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aM   a'M' 

where  a  is  the  speed  of  sound.   Cancelling  the  Mach  numbers 
and  substituting  N  (usually  RPM)  for  the  angular  speed  Q, 
one  obtains  the  relationship 

Nr  ^   N'r' 
a     a 

Since  the  compressor  performance   maps  usually  relate  the 
characteristics  of  a  specific  compressor  under  differing 
conditions,  the  radius  r  could  be  eliminated  and  the  similar 
rotational  speeds  expressed  as 


—  =   —r    ,    orN  =  N'  —7- 
a    a'  a' 


Since  a^  =  yRT,  the  ratio  of  sonic  speeds  reduces  to  the 
ratio  of  the  square  root  of  the  temperature  if  y  and  R  are 
the  same  for  both  compressors.   This  is  usually  tlie  case 
since  y  is  a  constant  for  air  within  the  temperature  range 
encountered  and  y  -  y'  is  fundamental  to  similarity.   The 
gas  constant  R  varies  with  molecular  weight,  however;  and 
if  the  fluid  medium  is  air  plus  foreign  gas  (steam  or  pro- 
ducts of  combustion  for  example) ,  then  the  similarity  rela- 
tionship must  retain  sonic  velocity  as  a  variable.   In  fact 
some  literature  (in  particular  older  literature)  uses  the 
corrected  rotational  speed  notation  of 

std        a 
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where  L  is  a  characteristic  length,  usually  the  tip  radius 
of  the  first  stage  if  comparing  different  sized  compressors, 

All  corrected  rotational  speeds  in  this  paper  are  the 
more  usual  N//9  since  foreign  gas  distortion  is  not  con- 
sidered.  References  2  and  3  analyze  the  effects  of  foreign 
gas  distortion  or  vorticity  and  how  it  affects  compressor 
performance . 
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APPENDIX  B 

Carta  and  Adamczyk  [Ref.  10]  define  a  stall  angle 
parameter  a  as  the  ratio  of  instantaneous  incidence  angle 
to  the  stall  incidence  angle  of  the  compressor  blade.   They 
further  define  the  incidence  angle  a  to  be  the  difference 
between  the  stagger  angle  X  and  the  inlet  flow  angle  S-,  . 
From  Figure  14 

a=3-,  -A,  or3-,  =  a  +  X 

and  3-1  differs  from  a  only  by  a  constant  X  that  is  a  func- 
tion of  blade  geometry.   It  will  not  be  inconsistent  then 
to  plot  a  distribution  of  inlet  flow  angle  B-,  vs  9*  instead 
of  o   vs  6*  as  in  Ref.  10. 

In  this  semi-quantitative  analysis  a  nominal  clean  in- 
let (undistorted)  value  of  45°  will  be  assumed  for  B, .   The 
changes  in  inlet  flow  angle  at  the  points  A,  B,  B',  C,  D 
and  E  of  Figure  17,  using  the  change  in  U   and  U^^  of  Figures 
17  and  21,  are  predicted  graphically  in  Figure  B-1.   The 
predicted  inlet  flow  angles  are  plotted  in  Figure  21  with 
the  stall  parameter  a  from  Ref.  10. 

Although  tlie  magnitude  of  many  of  the  variables  in  this 
semi-quantitative  analysis  are  somewhat  arbitrary,  the  di- 
rections in  which  they  were  applied  have  been  developed 
carefully  in  tlie  text  of  tliis  paper.   Deviations  from  the 
magnitude  used  slvould  give  tlie  s;uiie  general  shape  of  inlet 


67 


Point  A 


Point  B 


o 

•H 
4-> 

o 

G 

i-i 

o 
■p 
o 
^1 


6^  -  41° 


AU 


Point  B' 


=  0 


Point  C 


B-^  =  58° 


Point  D 


Point  E 


B^  =  40° 


=  0 


W.  for  nominal  clean  inlet  (45°) 
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flow  angle  curve  as  that 'of  Figure  21.   The  important  point 
is  that  the  general  trend  of  inlet  flow  angle  closely  parallels 
that  calculated  by  Carta  and  Adamczyk  and  reinforces  the  bent 
vortex  model  as  having  accounted  for  the  upstream  effect  of 
the  compressor  on  the  flow  field. 
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